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(HICs) from AGS to RHIC energies. The UrQMD hadron-string transport model as well as the 
CRAB analyzing program are adopted. Based on the cascade mode, in general, the calculations 
are satisfying and well in line with the experimental data although discrepancies are not negligible. 
Such as: I), the HBT time-related puzzle exists at all energies. II), at low AGS energies, the 
calculated volume as well as the mean free path of pion source at freeze-out are lower than data. 
It implies that a better description of interactions of particles at early stage of HICs is required. 
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1. Introduction 

For discovering the theoretically predicted quark gluon plasma (QGP) the heavy ions have 
been collided with nucleon-nucleon center-of-mass energies from less than y^s ~ 2.5 GeV (SIS/FAIR 
energy regime), 2.5 - 20 GeV (AGS/FAIR and SPS) up to 20 - 200 GeV (RHIC). Indeed, there 
are some signals - such as charmonium suppression, relative strangeness enhancement, etc. - of 
the (phase) transition to the deconfmed phase have been observed in heavy ion collisions (HICs) 
at SPS energies [III 0, [3L 0J|. Additional information about the space-time structure of the parti- 
cle emission source (the region of homogeneity) can be extracted by Femtoscopy ^ or namely 
Hanbury-Brown-Twiss interferometry (HBT) []6p. It is supposed that non-trivial structures in the 
excitation function of HBT quantities should be present at the energy threshold for the onset of 
QGP formation [Q]. Unfortunately, so far the excitation functions of the HBT quantities have not 
shown any obvious discontinuities within the large span of explored beam energies [||, |||9, 10, 11 
I|, 0, [B], m m 0> HI- 



A comprehensive theoretical investigation on the excitation function of the HBT parameters 
is thus highly required but still absent so far [|h. Recently, based on the UrQMD hadron-string 
transport model [^, 21 , |2|] and the CRAB [ ^3| , |4|] analyzing program, we investigated the beam 



energy, transverse momentum, system-size, centrality, and rapidity dependence of the HBT pa- 
rameters Ri, Ro, Rs (dubbed as HBT radii or Pratt radii), and the cross term Rol of pion source 
[25, 26, ^7|, 28]. In general, the calculations are satisfying and well in line with the experimental 



data although discrepancies are not negligible. Such as, I), the calculated Ri and R$ values for 
Au+Au collisions at low AGS energies are visibly smaller than the data if the default UrQMD ver- 
sion 2.2 (cascade mode) is adopted. II), the HBT-'puzzle' with respect to the 'duration time' of the 
pion source, is present at all energies. 

In this presentation, we show the beam energy and transverse momentum dependence of the 
pion HBT radii, as well as some other radius-related quantities such as the "duration-time" related 



quantity yjR 2 —R 2 s (and the Ro/Rs ratio), the freeze-out volume V/, and the mean free path A/ 
of pions at freeze-out. The standard UrQMD v2.2 in cascade mode is employed firstly to serve 
as a benchmark, then we find that the HBT time-related puzzle can be better understood with the 
consideration of a potential interaction. 

2. UrQMD transport model 

The UrQMD model is based on analogous principles as the Quantum Molecular Dynamics 
(QMD) [J||[3(|] and the RQMD [||] transport models. Similar to QMD, hadrons are represented 



by Gaussian wave packets in phase space, and the phase space of hadron i is propagated accord- 
ing to Hamilton's equation of motion: r, = and p, = — Here r and p are the coordinate 
and momentum of hadron i. The Hamiltonian H consists of the kinetic energy T and the effective 



two-body interaction potential energy U. In the standard version of UrQMD model [J20J, |2_1|], the 
potential energies include the two-body and three-body (which can be approximately written in 
the form of two-body interaction) Skyrme- (also called as the density dependent terms), Yukawa-, 
Coulomb-, and (optional) Pauli-terms as a base. Recently, in order to be more successfully ap- 
plied in the intermediate energy region (Eb < 2A GeV), more potentials are considered into the 
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model [32], those are, the density-dependent symmetry potential (essential for isospin-asymmetric 
reactions at intermediate and low energies) and the phenomenologically momentum-dependent in- 
teraction. We have found that the experimental pion and proton directed and elliptic flows from 
HICs with beam energies from ~ 100A MeV to 2A GeV can be well described with the poten- 



tials [J33j] . At higher energies, i.e., AGS/FAIR, the Yukawa-, Pauli-, and symmetry- potentials of 
baryons become negligible, while the Skyrme- and the momentum-dependent parts of potentials 



still influence the whole dynamical process of HICs. For example, in Ref. [|34|], with the help of a 
mean field from RQMD/S [ [35| ] and a Jet AA Microscopic Transportation Model (JAM) it has been 
found that the momentum dependence in the nuclear mean field is important for the understanding 
of the proton collective flows at AGS and SPS energies. 

In the latter part of the presentation, we will also show the importance of the mean field on the 
HBT time related puzzle. We adopt the soft equation of state (EoS) with momentum dependence 



(SM-EoS) which is same as that in Ref. Q34|]. Furthermore, as in [)34{] , the relativistic effects on the 
relative distance r,y = r, — r ; - and and the relative momentum p^- = p, — p ; (Lorentz transformation) 
employed in the two-body potentials are considered: 



ij 



J- ™2 



fijK*iJ HijJ 



Ej + Ej 



(2.1) 
(2.2) 



In Eqs. 2A and 12 the velocity-factor j3,- 7 and the corresponding /-factor of i and j particles are 
defined as 



P»- + Pj 

Ei+Ej 



, and Yij 



1 



■A 2 , 



(2.3) 



Similar to RQMD, the collision term of the UrQMD model treats 55 different baryon species 
(including nucleon-, delta- and hyperon- resonances with masses up to 2.25 GeV) and 32 different 
meson species, including (strange) meson resonances with masses up to 2.0 GeV as tabulated in 
the PDG as well as the corresponding anti -particles, i.e. full baryon/antibaryon symmetry 
is included. Isospin is explicitly treated as well. For hadronic continuum excitations a string 
model is used. Starting from the version 2.0, the PYTHIA is incorporated into UrQMD in order to 
investigate the jet production and fragmentation at RHIC energies [|22|]. 



3. Analyzing process 

To calculate the two-particle correlator, the CRAB program is based on the formula: 

c / k , = J^ 4 xi^ 2 g(^i,p 1 )g(x2,p 2 )|0(q,r)| 2 (3 ^ 

/ d 4 x i g {x i , p i ) / d 4 x 2 g (X2 , P2 ) 

Here g(x,p) is the probability for emitting a particle with momentum p from the space-time point 
x = (r,f). 0(q,r) is the relative two-particle wave function with r being their relative position, 
q = P2 — pi and k = (pj + p 2 )/2 are the relative momentum and the average momentum of the two 
particles. Due to the underlying quantum statistics, this correlator can be fitted approximately by 
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a Gaussian form. Using Pratt's three-dimensional convention (the LCMS system), the correlation 
function in Gaussian form reads 

C(qo,qs,q L ) = 1 + Xexp(-R 2 q 2 - R 2 q 2 - R 2 s q 2 s - 2R 2 OL q q L ). (3.2) 

Here qi and Rj are the components of the pair momentum difference q and the homogeneity length 
(HBT radii) in the i direction, respectively. The pre-factor A is the incoherence parameter and lies 
between (complete coherence) and 1 (complete incoherence) in realistic HICs. The term R 2 0L is 
called cross-term and vanishes at mid-rapidity for symmetric systems, while it deviates from zero 
at large rapidities [37, [27| ]. 



We compare our calculations of the Pratt parameters (at midrapidity) of the pion source with 
experimental data for the following centr al collisions of heavy nuclei: 

1. Au+Au at the AGS beam energies E b = 2, 4, 6, and 8A GeV (< 11% of the total cross 
section Oj), a rapidity cut |F cm | < 0.5 (Y cm = ^log( £ ""^p|| ), E cm and p^ are the energy and 
longitudinal momentum of the pion meson in the center-of-mass system) is employed. The 
experimental (E895) data are taken from [9]. 

2. Au+Au at the AGS beam energy 11. 6A GeV (the < 5% most central collisions), a rapidity 



cut \Y cm \ < 0.5 is employed. The experimental (E802) data are taken from [ 1C]. 



Pb+Pb at the SPS beam energies E b = 20, 30, 40, 80, and 160A GeV (< 7 '.2%a T of most 
central collisions), a pion-pair rapidity cut \Y nn \ < 0.5 (Y n7Z = ^log( ^|^^|||^^ ) is the pair 
rapidity with pion energies E\ and E% and longitudinal momenta pp and pn 2 in the center of 



mass system) is employed. The experimental (NA49) data are taken from Ql 1|, |i2Q. 



Pb+Au at the SPS beam energies E b = 40, 80, and 160A GeV (the < 5% most central col- 
lisions), the pion-pair rapidity cut Y n% = —0.25 ~ 0.25, —0.5 ~ 0, and — 1.0 ~ —0.5 are 



chosen. The experimental (CERES) data are taken from [13]. 



Au+Au at the RHIC energies = 30 (< \5%a T ), 62.4 (< \5%o T ), 130 (< 10%a r ), and 
200 GeV (< 5%o T ). Here a pseudo-rapidity cut \n cm \ < 0.5 (f] cm = ^logf^M, (p is the 
momentum of the pion) is employed. The experimental (PHOBOS, STAR, and PHENIX) 
data are taken from [0, 0, 00 |Hp. 



4. The HBT radii at AGS, SPS, and RHIC energies 

Fig. |] shows the kj (kj = (p ir +p 2r )/2) dependent HBT radii Ri (top), Ro (middle), and R$ 
(bottom) of the 7T~7z:~-pair for AGS energies 2, 4, 6, 8, and 11. 6A GeV (from left to right). The 
cascade mode is employed. We notice that the radii Rl and Rs are somewhat smaller than the 
data especially at lower beam energies, if a constant width for resonances (lines with diamonds) 
is used. Fig. |] also includes the results with the mass dependent lifetime of resonances (lines 
with triangles). With this treatment, the resonances with their small invariant masses decay later 
and hence the expanded fireball becomes larger as compared to the standard (mass independent) 
treatment. At large kj as well higher beam energies, this effect is reduced. It is interesting to see 
that the result with a mass-dependent treatment of resonance lifetimes can matches the data much 
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Figure 1: Transverse momentum kj 
dependence of the HBT-radii Rl, Ro, 
and R$ at AGS energies. The calcu- 
lated results with and without consid- 
ering the mass dependence of reso- 
nance lifetimes are demonstrated. Ex- 
perimental data are shown |l0| 
with scattered stars. 



better in the AGS energy region. It is also seen that this mass-dependence has almost no effect on 
the ratio between Ro and R$ values. 

In Figs. U and |3| we show the k T dependence of the HBT- radii R L (top plots), R (middle plots), 
and R$ (bottom plots) at SPS energies. In Fig. ||] the results are compared with preliminary NA49 
data [|l2|]. In Fig. || the results are compared with CERES data [13]. Furthermore, the 7i~ — %~ 
correlations are calculated in Fig. ^ while the two-charged-pion correlations (including two-7T~ 
and two-7T + mesons) are calculated in Fig. S. Firstly, it is very interesting to see that the present 
calculations can reproduce the kj -dependence of HBT radii Rl and R$ fairly well. Only at small 
kr, the calculated Rl and R$ values are seen up to 25% lower than data. While for the Ro values, 
the calculations are shown larger than both NA49 and CERES data especially at relatively large kj. 
By comparing the NA49 data with the CERES data for central Pb+Pb (o/o T < 7.2%) and central 
Pb+Au (a/ Or < 5%) collisions, one observes that the CERES Ro data are somewhat smaller than 
the recent published NA49 data [ |l2| ] especially at large kj and low beam energy Eb = 40A GeV 
although the recent NA49 data at large kj have already been driven down visibly when comparing 
with those preliminary data in the previous publication [11]. The origin of the difference between 
NA49 and CERES data is still not quite clear. At 160 A GeV (shown in Fig. ||) the HBT radii 
are also calculated by adopting zero formation time for strings (x s = fm/c). As a result, the 
kj -dependence of the HBT radii becomes steeper, and the values of R$ increase and approach the 
calculated Ro values. It should be noted that, although a shorter formation time is apt to explain 
the "HBT-puzzle", as well the elliptic flow, the absolute values of HBT-radii are not well in line 
with data. 

Fig. |] gives the kj dependence of the Pratt-radii Rl (left plots), Ro (middle plots), and R$ 
(right plots) at RHIC energies. Both the absolute values and the decrease of the Pratt-radii Rl and 
R$ with transverse momentum is reproduced by the present model calculations very well. Here, 
it is also seen that the calculated kj -dependence of R$ is somewhat flatter than that of Rl, which 
implies that flow effect on the kj -dependence of the Pratt-radii is still important. Besides the flow 
effect, the surface-like emission charactistic of microscopic models should play significantly role 
on HBT parameters as well because also other Cascade/Boltzmann model-calculations (see e.g., 
the RQMD [g g], the HRM [Q, and the AMPT @) can reproduce the k T dependence of Pratt 
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Figure 2: kj dependence of the Pratt-radii at SPS Figure 3: kj dependence of the Pratt-radii at SPS 
energies. Preliminary NA49 data are taken from energies. CERES data are taken from Jl3[], 
[0. At E h = 160A GeV the calculation results of 
the HBT-radii with a vanishing formation time for 
strings are also presented. 
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Figure 4: kj dependence of the Pratt- 
radii Ri, Ro, and R$ at RHIC ener- 
gies. Experimental data are shown 
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k T (MeV/c) 



radii (almost) equally well. The UrQMD calculations of Ri and R$ reproduce the experimental 
data well within the error bars, while the calculated Ro's are also larger than the experimental data 
— the Rp is about 25% too large. 



5. Excitation function of the Vf, the A/, and the ^R 2 — R 2 S (and Ro/Rs) of pions 

Fig. H (a) shows the excitation function of the pion source volume Vf at freeze-out, calculated 
as [|{J Vf = {2n)^RiR 2 s . The calculations with default cascade mode are shown at kj = 100 ± 
50 MeV (full line) and 200 ± 50 MeV/c (dotted line), respectively. The gray areas between the two 
lines are shown for better visibility. The data are at k T ~ 150MeV/c for AGS-E895 and SPS-NA49, 
at k T ~ 170MeV/c for reaction at ^/sm = 130 GeV, at k T ~ 200MeV/c for reactions at AGS- 



E802, SPS-CERES, and ^/s^ = 62.4, 200 GeV. Since the experimental data from NA49 [11,12] 
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Figure 5: (a): Exci- 
tation function of the 
Vf at transverse mo- 
menta between kj = 
lOOMeV and 200MeV 
(gray area), compared 
with data in this kj- 
region. (b): Excitation 
function of the X/ of pi- 
ons at freeze-out. (c): 
Excitation function of 
the duration-time related 
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Figure 6: Excitation function of the 
Ro/Rs ra tio at small kj. The data 
are indicated by solid stars. The 
dotted lines with open rectangles 
are results under cascade mode, the 
dashed lines with open diamonds 
represent the calculations with po- 
tentials of formed baryons. 



and from CERES [13] collaborations overlap at beam energies 40, 80, and 160A GeV, we show the 
calculations and data with respect to CERES energies separately as dashed-dotted lines and open 
symbols. Fig. || (a) shows clearly that the UrQMD cascade calculations do provide a reasonable 
freeze-out volume for the pion source at RHIC energies. At SPS energies, the agreement is fine 
with CERES data while it slightly underpredicts those of NA49. Towards even lower energies, the 
model underpredicts the measured freeze-out volume due to the omission of the strong interaction 
potential and other in-medium effects. E.g. at Eb = 2A GeV, the measured Vf is about 2 — 3 times 
larger than calculated value. As stated in Fig. [TJ, a mass-dependent lifetime of resonances accounts 
for an improvement of the HBT-radii at small kj and hence reproduce the data better. 
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The mean free path A/ of the pions at freeze-out is expressed as | |40| ] 

Yl . (5.1) 

NO NnOxn+NxOxx 

with the averaged pion-nucleon cross section On k = 72 mb and the averaged pion-pion cross section 
Gnu = 13 mb (Note that within the present model calculations these values are slightly energy 
dependent. However, here we have adopted the explicit numbers from Ref. p0| ] to compare to the 
results presented there). The nucleon and pion multiplicities Nn and N n are calculated as 

dN nuc i 

eons,pions 

~~dy ly "'"'"' * 



N N , n =y th W2n- ""T " lw (5-2) 



using the assumption of a thermal equilibrated system at freeze-out with a temperature Tf = 120 
MeV. Here, y t h is the estimated thermal homogeneity scale in rapidity at a certain kj and Tf, 
and is given by the expressions: y t f, = arctanh((/3,/,)), with (j3,/,) = yl + (j) 2 / (/) and (7) = 1 + 
1 /3 (Ki [mj /Tf) jKi(mj jTf) — 1) + Tfjmj. Here K n (z) is the modified Bessel function of order n 



and mj = \ m\ + kj. dN/dy\ } , mjcl is the rapidity density of pion (nucleons) at mid-rapidity. Recent 



calculations using the present UrQMD transport model [ |22| ] have shown that the calculated pion 
and nucleon yields are reasonably in agreement with data. 

Fig. H (b) shows the excitation function of Xf of pions at freeze-out. The experimental value 
for Xf at -y/SNN = 200GeV is obtained with the help of recent dN /dy data in plj], at all other 
energies the Xf data are taken from pC|]. It is seen that the theoretical Xf value increases gradually 
from ~ 0.5 to ~ 1 fm from AGS to highest RHIC energies with a weak dependence on kj. The 
experimental values of Xf are also between 0.5 — 1 fm. The observation (both experimentally and 
theoretically) of a nearly energy independent mean free path on the order of 0.7 fm at pion freeze- 
out is rather surprising. Physically it has been interpreted as a rather large opaqueness of the pion 
source at break-up [42, 



The HBT duration time "puzzle", i.e. the fact of the theoretical quantity ^R 2 —R 2 S being 
larger than extracted from the data, is present at all investigated energies (see Fig. || (c)): The 
calculated values of ^Jr 2 — R 2 s are about 3.5 ~ 5 fm while the measured ones are 1 .5 ~ 4 fm. Many 
efforts have been put forward over the last years to clarify this issue 

Here, we show in Fig. ^ the potential effect on the excitation function of the Rq/Rs ratio at the 



small kj. As seen in Fig. || (c), the Ro/Rs ratio, which is equivalent to the quantity t ~ y^o ~^S' 
is larger than the experimentally observed values at all investigated energies, if the cascade mode 
(dotted lines with open rectangles) is employed. When the SM-EoS is considered for the formed 
baryons (solid lines with open diamonds), the Ro/Rs ratio is seen obviously smaller than that 
with cascade mode and reproduces the (energy dependence of) data at AGS energies. At SPS and 
RHIC, however, the Ro/Rs ratio is still increasing monotonically with increasing beam energies 
and deviates from data again. At RHIC energies, the ratio approaches the one with cascade mode. 
It implies that the potential of formed baryons is increasingly losing its importance with increasing 
beam energies. At SPS and RHIC energies the deviation from data might be interpreted by the 
absence of the interactions of unformed particles from string fragmentation. Investigations are in 
progress. 
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6. Conclusion and Outlook 

To summarize, we show the transverse momentum and beam energy dependence of the HBT 



radii Rl, Ro> an d Rs> the quantity yR^—R^ (and the Ro/Rs ratio), the volume Vf, and the mean 
free path A/ of pions at freeze-out for heavy systems with energies from AGS to RHIC. In general, 
the model calculations with UrQMD v2.2 (cascade mode) are in line with the data over the whole 
inspected energy range. We also find a nearly constant mean free path for pions on the order of 
Xf = 0.7 fm which indicating a significant opaqueness of the source. 

Discrepancies especially in the lower AGS energy region are found and have to be resolved. 
The HBT duration-time related "puzzle" is present at almost all energies. The consideration of 
potentials for formed and unformed particles provides new insights into the origin of the time- 
related puzzle and the dynamics of HICs especially at early stage. 
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